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Abstract: We discuss two research areas dealing respectively with (1) a class 
of multivariate medians and (2) a symmetrization algorithm for probability 
measures. 



1. Introduction 

Geometric and stochastic ideas interact in a wide variety of ways over both theory 
and applications. For two connections that do not seem to have been mentioned 
before, we present here descriptive comments and suggestions for further investiga- 
tion. The first deals with multivariate medians, an active area of research in which 



Yehuda himself was interested (see [11| with Cun-Hui Zhang). Using the intrinsic 
volume functionals for convex bodies, we define a class of multivariate medians and 
show that among them, as special cases, are the L\ median and the Oja-median. 
The second question deals with a measure-theoretic generalization of the classic 
Steiner symmetrization technique for convex bodies. 



2. A class of multivariate medians 
2. 1 . Background 

In K d with its usual algebraic and metric structures, we consider the class JC d of 
convex bodies (compact, convex), which is closed under scaling 

XK = {Xx:xe K} 

and Minkowski addition 

K + L = {x + y : x g K, y G L}, 

A special class of convex bodies are sums of line segments, or zonotopes. The A- 
parallel body to K is K + XB, where Bd is the closed unit ball in M. d . Steiner's 
formula gives the volume of a typical parallel body: 

d 

vo\ d {K + XB d ) = V d -j(K) voljiBj) X j . 

3=0 
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The coefficients Vj(K), j = 0, 1, 2, . . . ,d are the so-called intrinsic volumes of K 
([HI). They arise in a variety of problems and can, for example, be defined quite 
differently: 

(2^^E^o\ J (Z bA K) 

Vj[K)= — j^m — ' 

where Z\j,d] , is a j x d matrix of independent, standard Gaussian random variables 
f[rH|). Some can be identified with familiar geometric functionals: 

V (K) = 1 

V\{K) = intrinsic width 



Vd-i(K) = 1/2 • surface area of K 
Vd(K) — d-dimensional volume of K 
Vj{K) :=0 for j > d. 

Here are some specific intrinsic volumes: 

V 1 ([0,1)) = 1 

n 

Vi([ai,6i] x [02,62] x • • ■ x a n ,b n ]) = /Xh - a { ) 

i=l 

and generally for j = 1,2, ... ,d: 

V,-([ai,&i] x [a 2 ,b 2 ] x ■•• x [a n ,b n ]) = ^ (6^ -a i:1 )(6 i2 - a i2 ) ■ ■ ■ (by - a^). 

il<i2<---<ij 



2.2. Vj -medians 

Suppose that a point sample x\,x 2 , ■ ■ ■ ,x n S R d is given, and a median is sought. 
We begin by creating, for every x £ M d , the zonotope 

Z(x) — x — x\ + x — x 2 + ■ ■ ■ + x — x n . 

Z(x) evidently aggregates, as a polytope, the discrepancies between x and the sam- 
ple points Xi. This can be quantified as follows: from our previous discussion, it is 
evident that the intrinsic volumes are generally measures of size. More specifically, 
Vj is a volume-like functional of homogeneity degree j. Consider the class of associ- 
ated variational problems: for each j = 1, 2, . . . , d, minimize Vj(Z(x)) over x 6 R d . 
The minimizing point (if it exists) we call the Vj -median of the sample. 

This point of view unifies two well-known medians that, on other grounds, appear 
to have little in common: 

• The V\ -median follows from minimizing ||x — Xi || and thus coincides with 
the well-known L\ median (an interesting treatment of computational issues 
and other properties appears in 
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• The V^-median, on the other hand, depends on minimization of vo\d{Z(x)) 
and thus coincides with the afhne-equivariant median of Oja [7J. 

Further work: 

• It would be of interest to investigate the Vj medians for intermediate values 
of j , 2 < j < d — 1. For such a j , one would seek to minimize 

^ ^ | det (x Xi-^ , x Xj 2 j . . . j x Xij ) | . 

l<ii <i2<-'-<ij <« 

• The Wills functional 

W{K) =1 + V!(K) + V 2 {K) + ■■■ + V d {K) 

is not homogeneous in K but on account of its other remarkable properties 
(e.g. [IH, [3], [13]) would likely be an interesting alternative size measure. 
The representation 

W{Z( X )) = f e -^ist 2 (y,Z(x)) dy 

could be useful. 

• One can focus on the inverse, or polar, body related to Z{x). One would then 
seek to maximize, for example, its volume, which is proportional to 

du. 



En , i-j 



• Finally, it would be of interest to relate these ideas to other instances of 
zonotopes' use in data analysis, e.g., @, (l3| . 



3. Steiner symmetrization as mass transport 

Steiner symmetrization is a geometric transformation of convex bodies that has 
been useful in a variety of problems, notably when an extremizing body is sought 
for a prescribed functional. Mass transport, on the other hand, embraces a number of 
issues that deal with shifts of mass (i.e., measure) from one location to another (see, 
for example, @, Q). Recently the question has arisen as to whether it is possible 
to generalize Steiner symmetrization along the lines of mass transport. Here we 
suggest a formulation and sketch some preliminary thoughts. 

Classical Steiner symmetrization is as follows: given a convex body K and a 
direction (= unit vector) u, locate each chord of K parallel to u and shift it along 
its line of inclusion so as to re-position its midpoint on the hyperplanc perpen- 
dicular to u. The aggregate of such shifted chords forms a new convex body, known 
as the Steiner symmetral of K. 

The original idea for this transformation was apparently due to L'Huiller in the 
1780s, but Steiner popularized it in his treatment of the isoperimetric inequality ([3(, 
[H). Variants of Steiner symmetrization can be found in [if, Q, The present one 
is motivated by re-casting the procedure in an equivalent form: regard each chord 
as bearing a uniformly distributed mass and shift it so as to re-position its center of 
mass on H^- . More generally one can think of a non-negative (and otherwise nice) 
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/ : R — > M , which defines a finite mass distribution on K and therefore on each 
line {u 1 - + tu, — oo < t < oo}, u 1 - 6 H^. On such a line, the center of mass is 

( [ /(u 1 + tu)dt] [ (u 1 - + tvjfiu 1 + tu)dt = u 1 - + mu. 

The shift then amounts to replacing t *—> f(u + tu) with t i— ► f(u + (t — m)u). 
The combined effect of such shifts (over w 1 ) amounts to a transport of mass. 

One can extend this even farther by regarding total mass as normalized to 
one and observing that the procedure has an equivalent formulation in terms of 
conditional expectations of random variables. This avoids issues of regularity for 
the density. Suppose then that X is a random vector in M. d for which EX exists 
( -E^ll-Xll < °°)- Then its Steiner symmetrization in the direction u e S^ 1 is 
defined to be the random vector 

(1) X u = X - E [uu'X | n„il] = X - uE [u'X | U U ±X] . 

Here II u x = I — uu 1 is (orthogonal) projection onto the subspace u . Thus, condi- 
tioned on H U ±X, the random vector X is shifted so that its conditional expectation 
lies inu 1 . 

It can be shown rather directly that this formulation extends classical Steiner 
symmetrization: 

Theorem 1. Suppose that X is uniformly distributed on a convex body K. Then 

(i) X u is uniformly distributed on the Steiner symmetral of K, and (ii) there is a 
sequence of symmetrizations that converge in distribution to uniform measure on 
the (centered) ball of the same volume as K. 

As we discuss below, general results appear to be difficult to prove. But the spe- 
cial case of symmetrization of a Gaussian measure already presents some interesting 
phenomena. We provide details for completeness. 

Theorem 2. Suppose that X has Gaussian distribution in R d with mean fi and 
(invertible) covariance matrix X, X ~ iV(/i, £). Then (i) any Steiner symmetriza- 
tion yields another Gaussian distribution, and (ii) there is a sequence of Steiner 
symmetrizations producing a limiting Gaussian distribution that is spherically sym- 
metric about the origin. 

Proof, (i) Direct by the linear nature of the transformation. 

(ii) There is no loss of generality in assuming [i = since, if not, a symmetrization 
with u = /VIImII yields this centering. For the first actual symmetrization, fix u and 
then as above 

X u = X-uE [u'X | IL U ±X}. 

For convenience, let II = H u ± . Well-known properties of normal random vari- 
ables provide that E [u'X \ IIX] is that linear functional of HX which minimizes 
E {u'X - w'nX) 2 . Setting VE [u'X - w'Ux] 2 = leads to 

nzn'w = nsw 

E TIw — T,u + cu for a constant c 
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For c, we apply vf to get = 1 + cm'E _1 u ==> c = — . Then 

n w = n 2 w = n [u + cE _1 u] = cim^u 

so that 



and 



E [u'X | HX] = -cu'YT^TLX 



X u = X - am'Tr 1 TLX 



(7 - cuu'E-^X. 



It follows that 

E„ = EX U X' U 

= E(I - cuu'T,- 1 U)XX'(I - cffirW) 

= (/ - cii«'s -1 n)E(/ - cinrW). 

Now suppose that (vi,Ai) and (v 2 , A 2 ) are two eigenpairs of E (||i>i|| = || 1^2 1 
A convenient choice for symmetrization is u = (wi + V2). Then 



11 = 1- -{v 1 +v 2 )(v' 1 +v' 2 ), 



and 



(2) 



1 



1 



1/1 1 



- = -u'E L u = --(vi + -y 2 )E~ 1 (wi + V2 ) = -- — + — 

Z Z V Al A 2 



Further, 



E -t n = + wa) ( -j-v[ + i« 2 ) ( 7 - >i + «a)K + v' 2 ) 



A 2 



1 



(vi + v 2 ) 



4 

1 ( 1 



1, 1, 1, 1, 1, 1 .' 

Ai A 2 Ai Ai A 2 A 2 

f («i-«2) + f («2-«i) 
Ai A 2 



1 

4 VAi A 2 
With lO, this gives 

cutt'E" 1 !! 



1 A 2 - Ai 

2 Ai +A 2 



(«1 +«2)K -V2). 
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Then 



1 A 2 - Ai 
" 2 Ai + A 2 

1 A 2 - Ai 
f 2 Ai + A 2 
/+ 1A 2 



(«1 

Ai 



2 Ai + A 2 
1 A 2 - Ai 



w 2 )(Aiu^ - A 2 w 2 ) 
Oi - v 2 ) (v[ + v 2 ) 



E II- 



1 A, 



2 Ai + A 2 



(vi - v 2 )(v' 1 +v' 2 ) 



2 Ai + A 2 
1 A 2 - Ai 



{vi + v 2 )(X 1 v' 1 - X 2 v' 2 )) + 



1 A 2 - Ai 

2 Ai + A 2 



(Ai«i - A 2 w 2 ) («i + v 2 ) 



S+ - 



2 Ai + A 2 
1 A 2 - Ai 



(wi + v 2 )(Xiv' 1 - A 2 w 2 )(ui - V2)(v[ + v' 2 ) 
[(vi + u 2 )(Aiui - \ 2 v' 2 ) + (Ai«i - A 2 w 2 ) (vi + u 2 ) 



2 Ai + A 2 

Xi)(vi+v 2 )(v[+v' 2 ) 



s + - 



1 A 2 - Ai 



2 Ai + A 2 



3Ai + A 2 , Ai + 3A 2 , 1 , 
2 o ~ 2^ ~ i^ 1 " 2 



Ignoring all but the first two terms of the spectral decomposition £ = XiViV^ 
\ 2 v 2 v' 2 H h \ d v d v' d , one gets 



Aiwi^ + A 2 w 2 w 2 + - 



1 A 2 



Ai 



2 Ai + A 2 



3Ai + A 2 , 

o ^1 



Ax + 3A 2 



-v 2 v 2 



A? + 6A!A 2 + \\ 



Xi){viv 2 + v 2 v[)} 
(v^ + v 2 v' 2 ) 



4(Ai + A 2 ) 
As a 2 x 2 array, this has eigenvalues 

Af+6A!A 2 + A2 



1 (A 2 
4 



Ai) 



Ai + A 2 



-{viv' 2 



v 2 v. 



± 



(A 2 -A!) 2 



4(Ai + A 2 ) 4(Ai + A 2 )' 



which simplify to 



1 



^(Ai + A 2 ), 



-(— — 

_2 Ui + A 2> 

It follows that E u has the same eigenvalues as £ except for Ai and A 2 , which 
are replaced by their arithmetic and harmonic means. Iterating this double-mean 
transformation leads to a common limit ([2|), i.e., Ai and A 2 can be replaced by a 
pair of eigenvalues arbitrarily close to one another. Since this is true for any pair 
of eigenvalues, the result is established. □ 

Remarks 

(1) The appearance of the two means in the proof is of historical interest, as the so- 
called "double-mean-iteration." See, for example, with a discussion of the work 
of Archimedes and Gauss among others. It would be of interest to know whether it 
appears here merely as a curiosity or whether there is a deeper theory connecting 
symmetrization of Gaussian measures and eigenvalue means. 

(2) In analogy with the classical theory, one would like to have results along the 
following lines: 
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Conjecture 1. For some sequence of Steiner symmetrizations [resp. for almost 
every sequence of Steiner symmetrizations (i.e. utilizing independent, random se- 
lections ofu, cf. jl4j )]- the induced sequence of symmetrized distributions is asymp- 
totically spherically symmetric. 

Conjecture 2. Any sequence of Steiner symmetrizations yields a (weakly) conver- 
gent sequence of probability measures. 

These do not seem easy to establish, but it seems plausible that results along 
the following lines may be helpful: 

Proposition 1. Steiner symmetrization is norm-reducing in mean-square: 

E\\X\\ 2 - E\\X U \\ 2 = E\\X - X u \\ 2 

= E[E(u'X\U u ±X)} 2 > 0. 

Proof. The following case is generic: d = 2, u = ( Q ) 

E\\X\\ 2 - E\\X U \\ 2 = E\\X ~ X u \\ 2 

= E[E{X 1 \X 2 )} 2 >Q. □ 

(3) Finally, it would be of interest to place the transformation X i — > X u in a 
martingale context. 
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